3.3.97 \(\int \frac {(a+\frac {b}{x})^n}{x^2 (c+d x)^2} \, dx\) [297]

3.3.97.1 Optimal result
3.3.97.2 Mathematica [A] (verified)
3.3.97.3 Rubi [A] (verified)
3.3.97.4 Maple [F]
3.3.97.5 Fricas [F]
3.3.97.6 Sympy [F]
3.3.97.7 Maxima [F]
3.3.97.8 Giac [F]
3.3.97.9 Mupad [F(-1)]

3.3.97.1 Optimal result

Integrand size = 20, antiderivative size = 133 \[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=-\frac {\left (a+\frac {b}{x}\right )^{1+n}}{b c^2 (1+n)}+\frac {d^2 \left (a+\frac {b}{x}\right )^{1+n}}{c^2 (a c-b d) \left (d+\frac {c}{x}\right )}-\frac {d (2 a c-b d (2+n)) \left (a+\frac {b}{x}\right )^{1+n} \operatorname {Hypergeometric2F1}\left (1,1+n,2+n,\frac {c \left (a+\frac {b}{x}\right )}{a c-b d}\right )}{c^2 (a c-b d)^2 (1+n)} \]

output
-(a+b/x)^(1+n)/b/c^2/(1+n)+d^2*(a+b/x)^(1+n)/c^2/(a*c-b*d)/(d+c/x)-d*(2*a* 
c-b*d*(2+n))*(a+b/x)^(1+n)*hypergeom([1, 1+n],[2+n],c*(a+b/x)/(a*c-b*d))/c 
^2/(a*c-b*d)^2/(1+n)
 
3.3.97.2 Mathematica [A] (verified)

Time = 0.26 (sec) , antiderivative size = 125, normalized size of antiderivative = 0.94 \[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=-\frac {\left (a+\frac {b}{x}\right )^n (b+a x) \left ((a c-b d) (a c (c+d x)-b d (c+d (2+n) x))+b d (2 a c-b d (2+n)) (c+d x) \operatorname {Hypergeometric2F1}\left (1,1+n,2+n,\frac {c \left (a+\frac {b}{x}\right )}{a c-b d}\right )\right )}{b c^2 (a c-b d)^2 (1+n) x (c+d x)} \]

input
Integrate[(a + b/x)^n/(x^2*(c + d*x)^2),x]
 
output
-(((a + b/x)^n*(b + a*x)*((a*c - b*d)*(a*c*(c + d*x) - b*d*(c + d*(2 + n)* 
x)) + b*d*(2*a*c - b*d*(2 + n))*(c + d*x)*Hypergeometric2F1[1, 1 + n, 2 + 
n, (c*(a + b/x))/(a*c - b*d)]))/(b*c^2*(a*c - b*d)^2*(1 + n)*x*(c + d*x)))
 
3.3.97.3 Rubi [A] (verified)

Time = 0.30 (sec) , antiderivative size = 150, normalized size of antiderivative = 1.13, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.300, Rules used = {1016, 948, 100, 25, 90, 78}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx\)

\(\Big \downarrow \) 1016

\(\displaystyle \int \frac {\left (a+\frac {b}{x}\right )^n}{x^4 \left (\frac {c}{x}+d\right )^2}dx\)

\(\Big \downarrow \) 948

\(\displaystyle -\int \frac {\left (a+\frac {b}{x}\right )^n}{\left (\frac {c}{x}+d\right )^2 x^2}d\frac {1}{x}\)

\(\Big \downarrow \) 100

\(\displaystyle \frac {d^2 \left (a+\frac {b}{x}\right )^{n+1}}{c^2 \left (\frac {c}{x}+d\right ) (a c-b d)}-\frac {\int -\frac {\left (a+\frac {b}{x}\right )^n \left (d (a c-b d (n+1))-\frac {c (a c-b d)}{x}\right )}{\frac {c}{x}+d}d\frac {1}{x}}{c^2 (a c-b d)}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {\int \frac {\left (a+\frac {b}{x}\right )^n \left (d (a c-b d (n+1))-\frac {c (a c-b d)}{x}\right )}{\frac {c}{x}+d}d\frac {1}{x}}{c^2 (a c-b d)}+\frac {d^2 \left (a+\frac {b}{x}\right )^{n+1}}{c^2 \left (\frac {c}{x}+d\right ) (a c-b d)}\)

\(\Big \downarrow \) 90

\(\displaystyle \frac {d (2 a c-b d (n+2)) \int \frac {\left (a+\frac {b}{x}\right )^n}{\frac {c}{x}+d}d\frac {1}{x}-\frac {(a c-b d) \left (a+\frac {b}{x}\right )^{n+1}}{b (n+1)}}{c^2 (a c-b d)}+\frac {d^2 \left (a+\frac {b}{x}\right )^{n+1}}{c^2 \left (\frac {c}{x}+d\right ) (a c-b d)}\)

\(\Big \downarrow \) 78

\(\displaystyle \frac {d^2 \left (a+\frac {b}{x}\right )^{n+1}}{c^2 \left (\frac {c}{x}+d\right ) (a c-b d)}+\frac {-\frac {d \left (a+\frac {b}{x}\right )^{n+1} (2 a c-b d (n+2)) \operatorname {Hypergeometric2F1}\left (1,n+1,n+2,\frac {c \left (a+\frac {b}{x}\right )}{a c-b d}\right )}{(n+1) (a c-b d)}-\frac {(a c-b d) \left (a+\frac {b}{x}\right )^{n+1}}{b (n+1)}}{c^2 (a c-b d)}\)

input
Int[(a + b/x)^n/(x^2*(c + d*x)^2),x]
 
output
(d^2*(a + b/x)^(1 + n))/(c^2*(a*c - b*d)*(d + c/x)) + (-(((a*c - b*d)*(a + 
 b/x)^(1 + n))/(b*(1 + n))) - (d*(2*a*c - b*d*(2 + n))*(a + b/x)^(1 + n)*H 
ypergeometric2F1[1, 1 + n, 2 + n, (c*(a + b/x))/(a*c - b*d)])/((a*c - b*d) 
*(1 + n)))/(c^2*(a*c - b*d))
 

3.3.97.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 78
Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(b 
*c - a*d)^n*((a + b*x)^(m + 1)/(b^(n + 1)*(m + 1)))*Hypergeometric2F1[-n, m 
 + 1, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m}, x] 
 &&  !IntegerQ[m] && IntegerQ[n]
 

rule 90
Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p 
_.), x_] :> Simp[b*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d*f*(n + p + 2))), 
 x] + Simp[(a*d*f*(n + p + 2) - b*(d*e*(n + 1) + c*f*(p + 1)))/(d*f*(n + p 
+ 2))   Int[(c + d*x)^n*(e + f*x)^p, x], x] /; FreeQ[{a, b, c, d, e, f, n, 
p}, x] && NeQ[n + p + 2, 0]
 

rule 100
Int[((a_.) + (b_.)*(x_))^2*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^( 
p_), x_] :> Simp[(b*c - a*d)^2*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d^2*(d 
*e - c*f)*(n + 1))), x] - Simp[1/(d^2*(d*e - c*f)*(n + 1))   Int[(c + d*x)^ 
(n + 1)*(e + f*x)^p*Simp[a^2*d^2*f*(n + p + 2) + b^2*c*(d*e*(n + 1) + c*f*( 
p + 1)) - 2*a*b*d*(d*e*(n + 1) + c*f*(p + 1)) - b^2*d*(d*e - c*f)*(n + 1)*x 
, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && (LtQ[n, -1] || (EqQ[n 
 + p + 3, 0] && NeQ[n, -1] && (SumSimplerQ[n, 1] ||  !SumSimplerQ[p, 1])))
 

rule 948
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_. 
), x_Symbol] :> Simp[1/n   Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*x)^ 
p*(c + d*x)^q, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p, q}, x] && NeQ 
[b*c - a*d, 0] && IntegerQ[Simplify[(m + 1)/n]]
 

rule 1016
Int[(x_)^(m_.)*((c_) + (d_.)*(x_)^(mn_.))^(q_.)*((a_) + (b_.)*(x_)^(n_.))^( 
p_.), x_Symbol] :> Int[x^(m - n*q)*(a + b*x^n)^p*(d + c*x^n)^q, x] /; FreeQ 
[{a, b, c, d, m, n, p}, x] && EqQ[mn, -n] && IntegerQ[q] && (PosQ[n] ||  !I 
ntegerQ[p])
 
3.3.97.4 Maple [F]

\[\int \frac {\left (a +\frac {b}{x}\right )^{n}}{x^{2} \left (d x +c \right )^{2}}d x\]

input
int((a+b/x)^n/x^2/(d*x+c)^2,x)
 
output
int((a+b/x)^n/x^2/(d*x+c)^2,x)
 
3.3.97.5 Fricas [F]

\[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=\int { \frac {{\left (a + \frac {b}{x}\right )}^{n}}{{\left (d x + c\right )}^{2} x^{2}} \,d x } \]

input
integrate((a+b/x)^n/x^2/(d*x+c)^2,x, algorithm="fricas")
 
output
integral(((a*x + b)/x)^n/(d^2*x^4 + 2*c*d*x^3 + c^2*x^2), x)
 
3.3.97.6 Sympy [F]

\[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=\int \frac {\left (a + \frac {b}{x}\right )^{n}}{x^{2} \left (c + d x\right )^{2}}\, dx \]

input
integrate((a+b/x)**n/x**2/(d*x+c)**2,x)
 
output
Integral((a + b/x)**n/(x**2*(c + d*x)**2), x)
 
3.3.97.7 Maxima [F]

\[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=\int { \frac {{\left (a + \frac {b}{x}\right )}^{n}}{{\left (d x + c\right )}^{2} x^{2}} \,d x } \]

input
integrate((a+b/x)^n/x^2/(d*x+c)^2,x, algorithm="maxima")
 
output
integrate((a + b/x)^n/((d*x + c)^2*x^2), x)
 
3.3.97.8 Giac [F]

\[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=\int { \frac {{\left (a + \frac {b}{x}\right )}^{n}}{{\left (d x + c\right )}^{2} x^{2}} \,d x } \]

input
integrate((a+b/x)^n/x^2/(d*x+c)^2,x, algorithm="giac")
 
output
integrate((a + b/x)^n/((d*x + c)^2*x^2), x)
 
3.3.97.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\left (a+\frac {b}{x}\right )^n}{x^2 (c+d x)^2} \, dx=\int \frac {{\left (a+\frac {b}{x}\right )}^n}{x^2\,{\left (c+d\,x\right )}^2} \,d x \]

input
int((a + b/x)^n/(x^2*(c + d*x)^2),x)
 
output
int((a + b/x)^n/(x^2*(c + d*x)^2), x)